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ABSTRACT 


The  potential  of  a  recently  developed  material,  pyrolytic  graphite, 
as  a  skin  structure  for  re-entry  vehicles  led  to  an  investigation  of  its 
thermal  and  mechanical  behavior  during  transient  heating  and  at  elevated 
temperatures. 

The  general  three-dimensional  equations  of  thermoelasticity  are 
developed  for  a  material  whose  properties  exhibit  the  same  anisotropic 
behavior  as  those  of  pyrolytic  graphite.  These  equations  are  then 
solved  exactly  for  a  two-dimensional  temperature  distribution  which 
closely  approximates  the  temperatures  measured  experimentally  in  pyrolytic 
graphite  plates. 

An  oxyacetylene  flame  apparatus  was  used  to  subject  l/8  inch  and 
l/4  inch  thick  pyrolytic  graphite  plates  to  transient  heating.  Measure¬ 
ments  were  made  to  determine  the  temperature  distribution  through  the 
thickness  and  along  the  length  of  the  plates.  An  experimental  technique 
was  developed  for  measuring  transient  strains  up  to  a  temperature  of 
1500*  F  using  a  Tuckerman  optical  strain  gage.  The  technique  was 
applied  to  measure  the  thermal  strain  variation  in  the  heated  pyro¬ 
lytic  graphite  plates. 

The  thermal  strains  calculated  from  the  analytical  solution  are 
compared  with  the  measured  strains  and  generally  found  to  agree  within 
experimental  error.  Calculations  of  the  stress  distributions  show  that 
maximum  stresses  of  -1400  psi  and  -3000  psi  were  produced  in  the  l/8 
inch  and  l/4  inch  thick  plates,  respectively.  Finally,  the  effect  of 
anisotropy  on  the  state  of  stress  is  evaluated. 
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elastic  constants 

E  elastic  modulus  psi 

6  shearing  modulus  psi 

h  thickness  of  plate  in 

L  length  of  plate  in 

l,m,n  direction  cosines 

T  temperature  *F 

u,v,w  displacements 


cartesian  space  coordinates  (i  ■  1,2,3) 
body-force  components  (i  -  1,2,3) 
surface-force  components  (i  ■  1,2,3) 

X  x^/L,  dimensionless  length  coordinate 

Z  x^/h,  dimensionless  thickness  coordinate 


Greek  Symbols 

a  coefficient  of  thermal  expansion 

£,7,1),  i  anisotropic  elastic  parameters  [defined  by  eqn. 

strain  components  (1  ■  1,2 . 6) 

v  Poisson's  ratio 

cp  Airy  stress  function  [defined  by  eqn.  (11)] 

0^  stress  components  (i  ■  1,2,. ..,6) 

T  time 


in/ in-#F 
(17)1 
in/  in 


psi 

sec 
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1.0  INTRODUCTION 

One  of  the  more  Important  problems  that  arise  during  re-entry  of  space 
vehicles  at  hypersonic  speeds  Is  caused  by  the  skin  temperatures  attained 
as  a  result  of  the  aerodynamic  heating.  Adiabatic  wall  temperatures  under 
these  conditions  can  exceed  temperature  limitations  of  materials  commonly 
used  in  the  manufacture  of  such  vehicles.  The  difficulty  Is  compounded  by 
the  fact  that  these  elevated  temperatures  are  often  accompanied  by  severe 
temperature  gradients.  Thermal  stresses  are  produced  by  these  large 
temperature  differences  and  may  be  an  Important  factor  In  determining 
structural  requirements. 

Although  thermal  protection  systems  for  present  space  vehicles  as  well 
as  those  in  the  near  future  can  be  designed  and  built  with  available  mate¬ 
rials,  the  performance  of  these  vehicles  (e.g.,  range)  will  be  limited  by 
weight  requirements.  The  possibilities  of  improving  performsnce  have  led 
»  to  the  developswnt  of  several  new  types  of  materials  possessing  promising 

potentialities.  One  such  material,  which  has  undergone  extensive  Investi¬ 
gation  during  the  past  two  years,  Is  pyrolytic  graphite.  This  material  Is 
one  of  the  strongest  structural  materials  known  for  use  at  very  high 
temperatures.  Its  highly  directional  thermal  properties  are  such  that  It 
Is  considered  to  be  a  thermal  Insulator  in  a  direction  normal  to  the 
surface  of  deposition  (basal  plane),  while  being  a  thermal  conductor  in 
a  direction  parallel  to  the  surface.  These  unusual  properties  of  pyro¬ 
lytic  graphite  may  be  regarded  as  a  combination  which  offers  the  advantages 
of  a  composite  heat  shield  and  skin  structure. 

Because  of  the  Interest  and  need  for  Information  regarding  the 
behavior  of  pyrolytic  graphite  during  transient  heating  and  at  elevated 
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temperatures,  a  research  program  was  undertaken  at  the  University  of 

California,  under  the  sponsorship  of  the  Lawrence  Radiation  Laboratory, 

to  investigate  in  particular  the  thermal  and  mechanical  behavior  of  a 

flat  plate  specimen  when  one  surface  was  subjected  to  parallel  flow 

heating  from  an  oxyacetylene  flame.  The  oxyacetylene  heating  apparatus 

was  developed  and  used  in  previous  work^’^*  and  found  to  provide 

2 

an  average  heating  rate  of  approximately  AO  Btu/ft  -sec.  Pyrolytic 
graphite  specimens  selected  for  the  tests  were  6"  x  A"  x  1/8"  plates 
and  6"  x  A"  x  l/A"  plates. 

The  experimental  phase  of  the  investigation  included  measurements 
of  the  temperature  distribution  in  the  plates  and  the  thermal  strain  at 
the  unheated  surface  of  the  plates  during  transient  heating.  An  earlier 
report*  describes  the  technique  employed  in  measuring  the  temperature 
distributions.  That  report  also  presents  an  analytical  solution  for 
predicting  the  two-dimensional  transient  temperature  variation  in  a 
thermally  orthotropic  plate  and  compares  this  solution  with  the  experi¬ 
mental  data  for  a  pyrolytic  graphite  plate.  In  the  present  report  a 
technique  is  described  for  measuring  transient  strains  up  to  a  tempera¬ 
ture  of  1500*  F  using  a  Tuckerman  optical  strain  gage.  The  technique 
was  developed  and  used  to  measure  the  thermal  strain  variation  of  the 
heated  pyrolytic  graphite  plates.  The  stresses  corresponding  to  the 
measured  strains  were  found  to  reach  maximum  values  of  -1A00  psl  and 
-3000  psi  in  the  1/8  inch  and  l/A  inch  thick  plates,  respectively. 

There  are  a  number  of  theoretical  analyses  dealing  with  the  problem 

2  3  A 

of  thermal  stresses  in  isotropic  plates.  ’  *  A  rather  general  solution 
* 

Note:  Numbers  in  superscript  denote  references  listed  in  Bibliography. 


to  this  problem,  which  extends  the  analysis  of  Reference  5,  was  recently  pre¬ 
sented  by  Gatewood.^  He  considered  the  three-dimensional  thermal  stresses 
In  a  moderately  thick  elastic  plate  In  which  the  temperature  Is  taken  as  a 
polynomial  In  the  thickness  variable  but  with  relatively  large,  although 
restricted,  gradients  with  respect  to  the  coordinates  of  the  plane  of  the 
plate.  However,  the  properties  of  pyrolytic  graphite  are  not  Isotropic, 
but  are  highly  directional  so  that  the  available  Isotropic  solutions  cannot, 
In  general,  be  applied.  A  thermoelastic  analysis  of  such  a  material  must 
proceed  from  the  generalized  Hooke's  law  for  anisotropic  solids. 

In  contrast  to  the  Isotropic  case,  there  appears  to  be  little  liter¬ 
ature  available  which  deals  with  thermal  stresses  In  anisotropic  plates. 

Pell^  has  considered  the  thermoelastic  behavior  of  an  anisotropic  plate, 
having  a  plane  of  elastic  symmetry  parallel  to  the  plane  of  the  plate, 
using  thin  plate  theory.  However,  he  was  unable  to  obtain  a  general 

Q 

solution  of  the  problem.  More  recently  Nowacki  presented  a  method  of 

24 

solution  utilizing  Galerkin's  functions  for  the  case  of  plane  stress 
in  an  orthotropic  plate.  Although  the  method  reduces  the  problem  to 
that  of  determining  three  functions,  the  mathematical  difficulties  appear 
to  remain  formidable.  Nowacki  also  obtained  a  solution  in  closed  form 
for  the  thermal  stresses  in  an  infinite  orthotropic  plate  with  a  constant 
temperature  prescribed  in  a  rectangular  region  in  the  Infinite  plane  and 
zero  temperature  prescribed  outside  this  region.  The  applicability  of 
this  solution  is,  of  course,  limited  by  the  Imposed  idealization. 

In  the  present  analysis,  thermal  stresses  are  considered  in  a  trans¬ 
versely  anisotropic  plate  which  not  only  exhibits  orthotropic  symmetry 
but  also  exhibits  synmetry  about  an  axis  normal  to  one  of  the  orthogonal 
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planes.  Such  symetry  Is  exhibited  by  pyrolytic  graphite.  As  usually  done 
in  thermoelastic  analyses,  the  effects  of  thermoelastic  dissipation  and 

9 

Inertia  are  ignored  so  that  uncoupled  quasi-static  thermoelastic  theory 
can  be  applied.  The  analysis  then  proceeds  from  the  general  three-dimen¬ 
sional  equations  of  thermoelastlcity  without  the  use  of  plane  stress  or 
thin  plate  theory.  The  formulation  of  the  solution  is  similar  to  that 
of  Boley  and  Tblin's^  for  isotropic  free  beams.  Although  the  solution 
presented  is  for  a  particular  two-dimensional  temperature  distribution, 
the  general  method  of  analysis  appears  applicable  to  plates  having  a  more 
general  temperature  distribution  and  a  higher  degree  of  anisotropy. 
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2.0  THERMAL  STRESS-STRAIN  DISTRIBUTION  IN  A  TRANSVERSELY  ANISOTROPIC  PLATE 
In  this  section  the  general  equations  of  theraoelasticity  are  developed 
for  a  material  whose  properties  exhibit  the  same  anisotropic  behavior  as 
those  of  pyrolytic  graphite.  These  equations  are  then  solved  in  three- 
dimensions  for  a  temperature  distribution  which  closely  approximates  the 
temperatures  measured  experimentally  in  a  pyrolytic  graphite  plate. 

2.1  General  Thermoelastic  Equations 

Solutions  to  general  three-dimensional  problems  of  thermoelastic¬ 
ity  must  satisfy  the  familiar  equations  of  equilibrium  and  the  conditions 
of  compatibility.  Additionally,  they  must  satisfy  the  prescribed  boundary 
conditions  at  the  surface  of  the  particular  body  considered.  These  equa¬ 
tions  in  rectangular  cartesian  coordinates  x^,  x^,  and  x^  are  given 
in  Reference  10: 

Equilibrium  equations 


3a,  3a,  3o, 

. _ 1  ,  _ O  ,  _ D 

5q  +  5q  +  3q 


5q  +  sq  +  5q 


a) 


d",  do.  do. 

sq  +  5q  +  5q  ■  x3 

where  X^,  X^,  are  the  body-force  components  in  the  x^,  Xj,  x^ 

directions,  respectively. 
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Compatibility  conditions 


a*22  s.j2  ' 


±t  +  £i  . 

dx32  dx22  ^5 


fv  !S  *S 

^x2  Sx32  ' 


h\  b  I  _  \  .^*5  \\ 

d*2^X3  \  "  ^1  ^2  ^*3  / 


(2) 


aV 


b  I  **U  d*5  .  ^6 


*«3 

2  dx.dx 
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a  /*4 

s: 


4  ^5  *6  \ 

-  +  5q-^  j 


Boundary  conditions 


OjJ  +  0^0  +  O^n  ■  Xx 


o2«  +  04n  +  afii  -  X2 


°3n  +  °5*  +  °4n  "  X3 


(3) 


where  X^,  X^,  ere  the  components  of  the  prescribed  surface  tractions 

In  the  x^,  x^i  directions  respectively,  and  where  J,m,n  are  the 
direction  cosines  of  the  outward-drawn  surface  normal. 

The  equations  presented  above  are  based  on  purely  mechanical  and 
geometrical  considerations  and  therefore  apply  generally  to  any  elastic 
solid,  whether  isotropic  or  anisotropic,  Isothermal  or  nonisothermal. 
Furthermore,  a  uniqueness  theorem^  assures  that  a  solution  which  satis¬ 
fies  (1),  (2),  and  (3)  constitutes  an  exact  and  unique  solution. 

For  the  most  general  case  of  an  anisotropic  elastic  solid,  the 

generalized  Hooke's  law  relation  contains  36  elastic  constants. **  However, 

12 

it  can  be  shown  by  arguments  based  on  thermodynamics  that  not  all  of 

these  are  Independent,  and  it  is  now  generally  accepted  that  there  are  21 

13 

Independent  elastic  constants  for  the  most  general  case  of  anisotropy. 

When  thermal  effects  are  taken  into  account,  the  Hooke's  law  relations 
contain  6  additional  Independent  coefficients  of  thermal  expansion.  The 
number  of  independent  elastic  and  thermal  constants  reduce  considerably 
when  a  solid  exhibits  structural  symmetry  with  respect  to  a  plane  or 
about  an  axis. 

In  particular,  pyrolytic  graphite  is  a  material  which  exhibits 
symmetry  with  respect  to  three  orthogonal  planes  and  also  about  an  axis 
which  is  normal  to  one  of  the  orthogonal  planes.  Love**  refers  to  a 
material  exhibiting  such  symetry  as  "transversely  isotropic."  However, 
the  slightly  more  descriptive  phrase,  "transversely  anisotropic,"  will 
be  used  here  in  referring  to  the  material.  It  is  shown  in  Appendix  I 
that  for  this  case  of  symmetry  the  reduced  generalized  Hooke's  law 
relation,  written  in  rectangular  cartesian  coordinates,  contains  5 
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Independent  eleetic  constants  and  2  independent  thermal  expansion  coeffic¬ 
ients. 

Consider  the  transversely  anisotropic  plate  of  Figure  1.  The 
rectangular  cartesian  coordinates  are  oriented  such  that  the  isometric 
planes  are  parallel  with  the  orthogonal  planes  of  symmetry.  The  x^x^- 
plane  is  chosen  to  lie  parallel  to  the  basal  plane  (referred  to  as  the 
"a"-directlon  in  pyrolytic  graphite).  Then  the  x^-axis  is  an  axis  of 
synuetry,  parallel  with  a  normal  to  the  basal  plane  ("c"-direction). 

In  this  case  the  generalized  Hooke's  law  is  written  as  (cf. Appendix  I) 


*1  “  all°l  +  a12°2  +  a13°3  +  alT 

€2  “  a12al  +  all°2  +  a13°3  +  alT 

*3  "  a13ffl  +  a13°2  +  a33a3  +  a3T 

*4  "  a44a4 


where  a66  -  2<an-a12>. 

These  equations  may  be  rewritten  in  terms  of  normally  measured 
physical  constants  as  follows:* 


* 

e.g.  | 


see  Reference  13  for  this  procedure. 
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c,  » 

-i-  (J,  - 

!34c  . 

• 

1 

Eu  1 

En  2 

E33  3 

e„  - 

.^2o1 

+  -L-0, 

-  , 

- 

2 

E11  1 

E11  2 

E33 

«„  m 

Vl3  a 

V13 

*  '  0- 

+  JL, 

3 

E33  1 

E33  2 

E33 

ff3  +  aiT 


a3  +  a3T 


(4) 


e.  ■  - —  o. 

4  G44  4 


*6  -  V12>  «6 


It  is  often  convenient  to  express  the  conditions  of  compatibility 
(2)  in  terms  of  the  stresses.  The  Hooke's  law  relation  given  by  equations 
(4)  permits  this  transformation.  For  example,  taking  the  first  of  equa¬ 
tions  (2): 

Appropriate  substitutions  from  equations  (4)  yield: 


jL.v  _£L+jl(jL.v  jLL 
eiiU22  12  V  J  1  En  l  V  Vl2ax22r2 


(5) 
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Differentiating  the  first  and  second  of  the  equilibrium  equations  (1), 
taking  the  body-force  components  to  be  zero,  gives 


z\  A5  *6  °\ 

\2'  *1^3  ’  c*ldx2  ’  '  5xj2  ' 

Adding  and  using  the  third  of  equations  (1)  results  in 

d2o&  d2^  b2a2  ^  I  bo^  bo5  ^  d2*^  d2o2  d2o3 

2  *F5’  ‘  S? '  5^ '  ’  i^  +  ^7 

Finally,  substitution  of  this  result  into  (5)  leads  to 


By  following  similar  procedures  on  the  remaining  compatibility  conditions 
(2),  the  following  relations  are  obtained: 
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2  2V 

sqs;  (Ii +  V +  L  v  ‘  *7 


13 
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a2°3 

^x75x 
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(6d) 


j  ^  F  1  ^13  1  ^  °3 

'IT5!5! <<l +  °2>  +  [g^  -  J  sqs; 

+ fe  ( + i?)+  ^  <14Vi2>  *7] 


(6c) 


"5 +  “1  sqs^  ■ 0 


,G44  E33  .  ^*1^*2  ^  1  E33  dxldx2 


{e!*(<fc1*+^Ii)+*!i <l+Vi2>  42] 


(6f) 


d2T 

°6  +  *3  5^  “  ° 
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The  six  relations  given  by  equations  (6)  represent  the  general 
stress  formulation  of  the  compatibility  conditions  for  a  transversely 
anisotropic  material  with  constant  properties. 

2.2  Solution  for  a  Two-Dimensional  Temperature  Distribution 

Consider  again  the  transversely  anisotropic  plate  of  Figure  1. 
Let  the  plate  be  subjected  in  its  interior  to  a  temperature  distribution 
given  by 


T(xlfx3)  -  axL  +  bxx2  +  cxj,2x32  +  TH(x3)  (7) 

x3)  is  an  arbitrary  function 
of  x3<  The  particular  form  of  (7)  was  chosen  because  the  experimentally 
measured  temperature  distributions  were  found  to  be  closely  approximated 
by  such  an  expression  (cf.  Section  4.1).  The  plate  is  supposed  free  of 
external  loads  and  body  forces. 

The  stress  distribution  due  to  T^(x3)  in  a  transversely  aniso¬ 
tropic  plate  is  identical  to  that  in  an  isotropic  plate  with  the  same 
temperature  distribution.  This  is  attributed  to  the  isotropic  behavior 
of  the  transversely  anisotropic  plate  in  directions  parallel  to  its  plane. 
The  solution  for  this  case  is  given  by  Timoshenko^  who  arrived  at  the 
solution  by  physical  arguments  involving  the  application  and  removal  of 
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boundary  forces.  The  same  solution  is  derived  by  Boley  and  Weiner  , 
who  used  a  more  rigorous  semi-inverse  approach.  The  solution  leads  to 
the  following  expressions  for  the  stresses  in  the  plate: 


where  a,  b,  and  c  are  constants  and  T^( 


13. 


°3  ■  °4  ’  *5  ■  "6  ■  0 


if  1  12  (Xj-  2*1)  "| 

yi"^-^L'aiEiiT+^T+ — 


(8) 


where  h  is  Che  plate  thickness  and 


h  h  - 

NT  -  /  Q!1E11T  dx3  ;  Hp  -  /  Qt1E11T(x3-  jh)  dx3 


Obviously,  the  solution  is  not  correct  at  the  edges  of  the  plate  since 
the  stresses  do  not  vanish  there.  However,  by  Saint  Venant's  principle 
the  stresses  at  distances  away  from  the  edges  which  are  larger  than  the 
thickness  of  the  plate  will  be  correct  to  a  high  order  of  magnitude. 
Note  that  the  stresses  given  by  equations  (8)  take  into  account  the 
variation  of  properties  with  temperature. 

The  strain  components  corresponding  to  the  stresses  are  ob¬ 
tained  by  substitution  of  equations  (8)  into  the  Hooke's  law  relations 
given  by  equations  (4).  The  results  are 


« 


1 


1  „  ,  ll<*3- 

h  NT  +  73 


■ 


V13  1 

= — TTT? — r  CK.E..T  -  4  N_  - 
£33X^12  1  11  h  T 


12  (X, 


5*> 


+  a3T 


(9) 


«/  ■  «. 


-•e- 
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Because  of  Che  linearity  of  the  thermoelastic  equations,  the 
solution  given  by  equations  (8)  and  (9)  can  be  simply  added  to  the 
solution  which  follows. 

The  formulation  of  the  stresses  due  to  the  remaining  terms 

in  (7)  remains  to  be  developed.  The  analysis  to  follow  is  an  extension 
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of  the  exact  theory  of  Boley  and  Tolins  for  isotropic  free  beams. 

The  solution  will  be  obtained  by  the  semi-inverse  method  in 
which  the  posteriori  assumptions  are  made  that 


*1  *2  \  \  „ 
5^  -  ■  0 


(10) 


As  the  analysis  proceeds,  it  will  be  shown  that  all  the  general  equations 
of  thermoelasticity  developed  in  Section  2.1  can  be  satisfied  on  this 
basis. 

An  Airy  stress  function,  assumed  Independent  of  Xg,  is  now 
introduced  such  that 


Then  the  first  and  third  of  the  equilibrium  equations  (1)  are  satisfied. 
Substitution  of  (11)  into  the  compatibility  condition  (6e)  gives 


„  En  a2  /a2®  \ 


+  V 


12 


-  Vn 


d2T 


where  T  ■  T(Xj,x3) 
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Integration  with  respect  to  x^  and  yields 

°2  “  V13  ^2  4  V12  "  aiEllT  +  8l(xl)  +  g2(x3*  (12) 

Substitution  of  (12)  into  compatibility  conditions  (6a),  (6b),  and  (6c) 
leads  to  the  following  results: 


d*2 


dx  2 


>  (13) 


&.+ .-2  .  !u/f33-.v  .  v  v  L.  sjg_ 

3x 34  1-V122  E33  V  2G44  13  12  13  /  cbc^2 


-1--  hlL.y  2  hi  )£sl. 

1-V122E33\  13  E33/Sx14 


(14) 


<2.  En  ;>2_  a,  E„  /  E. .  a,  \n2_ 

.  .  -J-1I  (1+w  )  ix  .  XJ1  »  1>  +  X  XL 

1-V12  S  1-V122\13E33  Vdxj2 


Equations  (13)  give 


81  "  C1X1  +  c3* 


82  "  C2X3  +  C3" 
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Or 


gj  +  g2  *  +  c2x3  +  c, 


Thus  (12)  may  be  written  as 


E  2  2 

’2  '  ”13  *£  '  **  +  Vl2  '  W  ♦  Vl  +  «A  *  C3 


(15) 


1  “"3 

The  expression  given  by  (14)  will  be  abbreviated  in  the  following  forms 


^se_  +  p  -..fa  +  y  <£se_  .  <£x_  _5  £i_ 


^X^  ^Xj^dx-j2  dXj 


Sx, 


(16) 


where: 


P  ■ 


1-V 
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hil  hi.  _v  .  v  v 

E33\  2G44  13  1213 


-JL_  .  Ill  [  l-v  2  !n 

1-V122  E33  l  13  E33 


al  E11 
1  -  V 
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■  1  /  v  !il  +  ^3  ] 

i-v122  l  Vl3  e33  «1  / 


>  (17) 


To  solve  (16),  it  is  assumed  that 


q>  -  <p.  +  <p  +  .  .  .  -  2  <p 

i-1  1 


(18) 
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where  the  <p^'s  ere  chosen  as  follows: 


d4<Px 


52t 

ax2 


(19a) 


a4<p2 


5^2 


(19b) 


a*J 


D  *S 

P  ax1zSic32  bxf 


(19c) 


^i  .  n  ^*1-2 

dx^dxj2  dx^ 


(1  ■  4,5,6, • . . ) 


(19d) 


By  substitution,  it  is  easily  shown  that  equations  (19)  satisfy  the 
differential  equation  (16).  It  follows  from  (11)  and  (18)  that  the 
stresses  will  be  given  by 


al  " 


i  +  th.+ 


dx„ 


dx. 


°5  ' 


A  .  .  *\  .  sS 

Sx1<3x3  "  dx1^x3 


see 


(20) 
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a2q>.  a2<p 

- ?  +  - 7  +  * 

bx*  SXj2 
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The  following  boundary  conditions  hold  for  a  free  plate: 


o,  ■  — ■  0  at  x.  ■  0.  x,  ■  L 

1  d*32  11 


5  - 


a2 


^3 


§-  ■  0  at  x3  -  0,  x3  -  h 


d3  ■  —  2  "  0  »t  *3  ■  0*  X3  ■  h 

**1 


(21a) 


(21b) 

(21c) 


However,  It  Is  often  not  possible  to  satisfy  exactly  condition  (21a).  In 
such  cases  the  following  integral  boundary  conditions  will  be  substituted: 


h 

;  «! 

o 


h 

/ 


o 


h 

/  x3  dx3 
o 


h 

/ 


o 


0 


(21a*) 

(21aH) 


These  boundary  conditions  will  lead  to  the  correct  solution  within  the 
approximation  of  Saint  Venant's  principle.  Each  of  the  cp^'a  will  be 
required  to  satisfy  the  conditions  given  by  equations  (21). 

The  particular  temperature  distributions  to  be  considered  here 

x3)  so  that 

2  2  2 

T(xlfx3)  •  aXj  +  bxj  +  cxj  x3 


are  the  terms  in  (7)  other  than  T^( 


Substitution  into  (19a)  gives 


2  cqxj 


2 
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Integration  yields 


a*  2 


-  +  x3  SjCxj^)  +  S2(x3) 


where  Sj  end  Sg  are  undetermined  functions  of  x^.  It  is  apparent 
from  Inspection  of  this  result  that  the  exact  boundary  conditions  (21a) 
cannot  be  satisfied.  Therefore  (21a')  and  (21a")  are  applied 


i*  a  b  r  2  2  -i 

/  — 1  dx3  -  /  L  -  cHxt  x3  +  x3  S1(x1)  +  S2(x1)J  dx. 


3  2 

-  j-  clJXj2  +  S1  +  hS2  -  0 


^  ^  ^  r  2  2  2  i 

-  *3  i*1  6X3  '  -  L  ‘  °'*1  ‘3  +  *9  +  ‘3  S2<*3>  J  d,3 


J>  j  k3  k2 
—  cHXj  +  3  Sj  +  f"  S2  -  0 


Upon  solving  simultaneously 


Sj  -  cT|hx^ 


c  1  ,,2  2 

S2  -  -  g-  cljh  xx 


Thus 


(22) 
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Now  Integration  of  (22)  leads  to 


2  2  T  1  2  1  x3  i  x3  1 

*1  -  cijh  xx  -  12  *3  +  6  "h"  "  12  T?  J  +  *3  S3(xl)  +  W 


Applying  (21c)  results  In 


d2S, 


■f  -  0 


d2S. 


dx. 


-  0 


Hence 


Then 


and 


S3  ■  C4X!  +  e5  »  SA  -  c^x,  +  c 


4  6  1  T  7 


+  c4x3Xj  +  c5x3  +  c6»j  +  c7 


So  that 


d<pi  2  r  i  iv  ivi 

dxjdxj  *  2ct,h  X1  ’  6  *3  +  2  h  3  h2  +  C< 


Conditions  (21b)  denands  that 


Thus 


c.  ■  0 
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- i x3 ♦  |  ( 


(23) 


(24) 


21. 


From  (19b)  and  (23) 


^2  -  -  2b£  -  2c£xj2  -  2cMh*[- 


Upon  integrating 


d2cp 


_ 2  .  ^  2  1  »  4  „  2  1 

r~2  "  '  b5x3  '  6°Z  *3  "  2cpT|h  '  12 


x  2  +  i!l 

3  6  h 


4n 

_1  ^3_1 

12  h2J 


+  x3T1(x1)  +  T2(Xj) 


Applying  boundary  conditions  (21a1)  and  (21a")  to  this  result  and  solving 
for  and  gives 

Tl  -  b^h  +  ~  c£h3 

*  1  .  2  1  ...4  1  k.  4 

T2  -  -  g  b$h  -  Ygo  cPnh  -  30  c£h 

Thus 


Integration  of  (25)  yields 
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+  x3T3(x1)  +  T4(Xj) 
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Boundary  condition  (21c)  leads  to  the  conclusion  that 


d2T„ 


dx 


t  "  °  * 


d2TA 
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It  follows  that 


T3  "  Dl*l  +  D2 


I4  -  Vl+D4 


Thus 


d2q>2 


(26) 


and 


<P2  -  f(*3)  +  D1*1x3  +  D2x3  +  D3xx  +  D4 


a2<pz 

Sx^Sx3 


Hence 
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Condition  (21b)  demands  that  ■  0.  Therefore 


^2«p2 


3 


0 


(27) 


Now  from  (20c)  and  differentiation  of  equations  (23)  and  (25) 


n4 
°  *3 


Integrating, 


aV 

- j  -  x3U1(x1)  +  U2(Xj) 

dx- 


In  this  case  the  exact  boundary  conditions  (21a)  can  be  satisfied. 
Then 

Uj  -  U2  -  o 


So 


-  0 


(28) 


It  follows  that 


sqs;  ’  0 


(29) 


Thus  the  series  in  the  cp^'e  may  be  truncated  after  two  terms  since  the 
second  derivatives  of  the  remaining  terms  are  identically  zero.  The 
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expressions  for  end  cp2  constitute  the  solution  since  they  satisfy 
the  differential  equation  (16)  and  the  boundary  conditions  (21).  Now 
recall  that  0 2  is  given  by  (15)  as 


°2 

a 

*13  e33  +  Vl2  if  '  °lEuI  +  *  C2'3 

To  determine 

ci 

,  c^t  00^  C3  the  following  boundary  conditions 
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0 

x302dx3dx1  -  0 

The  following  result  is  then  obtained: 


-  VU  <*3 +  iio  "|h6> +  vi2°i ♦  “iEubl2  [-i  +  r-(r 


(30) 


Complete  expressions  for  0^,  o^,  0^»  and  0^  are  now 
evallable  by  appropriate  substitutions  in  equations  (20)  and  (30). 
Note  that  these  stresses  were  found  independent  of  0^  and  0^. 
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Now  and  will  be  required  to  satisfy  the  stress-free  boundary 
conditions  exactly.  Then  it  necessarily  follows  from  assumptions  (10) 
that  ' 

«4  '  C6  •  0 


The  remaining  compatibility  equations  (6d)  and  (6f)  and  the  second  of 
the  equilibrium  equations  (1)  are  identically  satisfied  by  this  choice. 

For  convenience  the  expressions  for  the  stresses  are  summarised 

below: 
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The  strains  corresponding  to  these  stresses  are  easily  obtained 
by  substitution  o£  equations  (31)  into  the  Hooke's  law  relations  (4).  In 
particular,  the  strain,  6^,  at  the  unheated  surface  (Z  ■  0),  which 
will  be  of  interest  in  comparing  with  experimental  data,  is  given  by 

-  ^-[1  c,|h2l'2x2  +  6  b£l,Z  +  30  c£h4  +  I5o  cWh4] 

■  TSo  -nrf  ^  +  vi2“lbI-2  [  f  -  *]  <32> 

+  V12011ch2L2  [|  -  -J  x]  +  C^aLX 

+  a1bL2x2(l-H»12) 

2.3  Calculated  Results 

The  solution  presented  in  the  previous  section  for  a  transversely 
anisotropic  plate  was  used  to  calculate  the  thermal  stresses  and  strains 
in  the  pyrolytic  graphite  plates  of  the  two  sizes  used  in  the  experimental 
program  (cf.  Section  3.2).  The  strain  at  the  unheated  surface  was  calcu¬ 
lated  from  equations  (9)  and  (32),  and  the  stresses  from  equations  (8)  and 
(31). 

It  should  be  recalled  that  equations  (31)  and  (32)  were  based 
on  constant  properties,  while  equations  (8)  and  (9)  accounted  for  variable 
properties.  Stresses  calculated  from  equation  (8)  were  found  to  be  sig¬ 
nificantly  affected  by  the  temperature  variation  of  the  thermal  expansion 
coefficient,  Of^,  in  the  "a"-direction  of  pyrolytic  graphite  (see  Figure 
24).  When  a  constant  average  value  of  was  used  in  equation  (8), 
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Che  resulting  stresses  were  found  to  differ  by  as  much  as  50  percent  from 
the  calculated  stresses  based  on  a  variable  a^.  A  similar  Influence  was 
found  on  the  strains  calculated  from  equation  (9).  These  calculations 
emphasised  the  Importance  of  accounting  for  property  variations  with 
temperature.  However,  although  equations  (31)  and  (32)  did  not  account 
for  variations  in  a^,  they  added  only  minor  contributions  to  the  stresses 
and  strains  so  that  the  errors  induced  by  the  use  of  a  constant  average 
value  of  Ctj  in  these  equations  were  felt  to  be  relatively  unimportant. 

Details  of  the  calculation  procedures  are  presented  In 
Appendix  II.  The  results  will  be  presented  and  discussed  In  Section  4.0. 
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3.0  EXPERIMENTAL  PROGRAM 

The  objectives  of  the  experimental  program  were  to  determine  the 
temperature  and  thermal  strain  distributions  in  pyrolytic  graphite  plates 
subjected  to  transient  heating  from  a  high- temperature  oxyacetylene  flame. 
Descriptions  of  the  material  tested,  the  heating  apparatus,  the  experi¬ 
mental  techniques,  and  the  test  procedures  will  be  presented  in  the 
following  sections. 

3.1  Material  Tested 

Pyrolytic  graphite  is  produced  by  the  deposition  of  a  gaseous 

hydrocarbon  onto  a  mandrel  of  commercial  graphite  which  is  maintained  at 

a  temperature  above  3600*  F.  The  properties  of  the  resulting  material 

can  be  varied  considerably  by  variations  in  the  processing  parameters. 

There  is  also  some  evidence  that  property  variations  may  occur  within 

an  individual  sample  depending  on  the  amount  of  residual  stress  present 
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in  the  material.  Recent  reports  indicate  that  it  may  be  possible  to 
classify  pyrolytic  graphite  into  as  many  as  16  different  types  based  on 
the  micro-structure. 

The  pyrolytic  graphite  specimens  selected  for  the  tests  were 
6"  x  4"  x  1/8"  plates  and  6"  x  4"  x  l/4"  plates.  Four  flat  plate  speci¬ 
mens  of  each  thickness  were  used  in  the  experiment. 

3.2  Heating  Apparatus 

The  heating  apparatus  used  for  the  experimental  phase  of  this 
investigation  was  developed  and  used  previously  for  the  determination 
of  the  behavior  of  ablating  materials  at  low  heat  rates.  A  detailed 
description  of  this  apparatus  and  its  heating  characteristics  are 
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given  in  References  16,  17,  and  18.  Only  the  important  features  will  be 
discussed  here. 

The  primary  elements  of  the  heating  apparatus  are  the  oxygen 
and  acetylene  fuel  sources,  the  fuel  lines  and  regulators,  and  multiple- 
nozzle  flame  head,  the  test  section,  and  the  water-cooled  exhaust  duct. 

The  general  arrangement  in  the  test  area  can  be  seen  in  Figure  2.  The 
apparatus  provided  an  average  heating  rate  of  approximately  40  Btu/ft  -sec 
to  one  surface  of  a  flat  plate  specimen  mounted  parallel  to  the  flame. 

The  velocity  was  300  ft/sec,  and  the  gas  temperature  approximately  5000*  F 

at  atmospheric  pressure.  Heating  rates  were  reproduced  to  within  about 

2 

1  2  Btu/ft  -sec  from  teat  to  test. 

For  the  purposes  of  the  present  investigation,  it  was  necessary 
to  modify  the  test  section  in  such  a  way  that  a  flat  plate  specimen  could 
be  mounted  at  the  top  of  the  section.  This  was  accomplished  by  surrounding 
the  sides  and  bottom  of  the  test  section  with  refractory  bricks  which  were 
lined  with  thin  plates  of  ordinary  connercial  graphite  to  alleviate  the 
flame  damage.  The  plate  was  then  placed  on  the  bricks  which  formed  the 
sides  of  the  section.  The  open  sections  in  front  of  the  leading  edge 
and  behind  the  trailing  edge  of  the  plate  were  covered  with  refractory 
bricks  so  that  the  oxyacetylene  flame  was  completely  enclosed.  Details 
of  the  test  section  assembly  are  shown  in  Figure  3.  With  the  plate  in 
position,  the  leading  edge  was  about  1/8"  above  and  l/2"  downstream  of 
the  top  row  of  nozzles  and  a  3  inch  x  6  inch  area  of  one  surface  of  the 
plate  was  exposed  to  heating  by  the  flame.  Flame  leakage  was  prevented 
by  filling  openings  with  Sauereisen  Ho.  7  and  by  maintaining  a  smooth 
fit  between  the  plate  and  the  bricks. 
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3.3  Technique  and  Test  Procedure  for  Temperature  Measurements 

Since  the  technique  and  test  procedure  for  measuring  the  transient 
temperature  distribution  in  a  pyrolytic  graphite  plate  have  already  been 
described  in  Reference  1,  only  a  brief  description  will  be  presented  here. 

The  basic  technique  was  to  insert  thermocouples  into  small  (1/8" 
diameter)  flat-bottomed  holes  milled  in  the  plate.  Starting  near  the 
unheated  surface,  these  holes  were  made  progressively  deeper  so  that  the 
temperature  variation  through  the  thickness  of  the  plate  could  be  measured. 
The  effect  of  the  holes  on  the  temperature  distribution  was  believed  to  be 
small  because  of  the  high  thermal  conductivity  of  pyrolytic  graphite  in 
its  basal  planes.  The  thermocouples  were  threaded  into  0.10  inch  diameter 
ceramic  tubes  and  spring-loaded  into  the  holes  with  the  device  shown  in 
Figure  A.  Number  32  gauge  chromel-alumel  thermocouples  were  used  and  a 
continuous  record  of  their  output  obtained  on  a  twelve -channel  Mlnneapolis- 
Honeywell  recording  oscillograph. 

Temperature  distributions  were  measured  in  both  the  l/8  inch 
and  l/4  inch  thick  plates.  Temperatures  at  the  unheated  surface  of  a 
plate  were  measured  first.  In  successive  tests  temperatures  were  measured 
in  holes  located  along  the  center-line  of  the  plate  parallel  to  the  flow 
at  distances  of  1/2,  1-1/8,  1-7/8,  3-3/8,  and  5-3/8  inches  from  the  leading 
edge  by  the  technique  described  above.  The  l/8  inch  thick  plate  was  heated 
for  a  period  of  approximately  30  seconds  in  each  test;  while  the  l/4  inch 
thick  plate  was  heated  for  approximately  70  seconds  in  each  test. 
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3.4  Development  of  Technique  for  Transient  High  Temperature 

Strain  Measurements 

During  the  early  phasee  of  this  investigation,  a  study  was  made 
of  the  available  experimental  techniques  for  measuring  transient  strains 
at  high  temperatures.  Of  the  various  techniques  studied,  the  electric- 
resistance  foil  gage  and  the  Tuckerman  optical  strain  gage  appeared  the 
most  promising.  Therefore,  preliminary  tests  were  conducted  using  both 
types  of  strain  gage. 

Foil  gages  were  bonded  to  pyrolytic  graphite  using  both  Trans- 
Sonlcs,  Inc.  TyP®  64-C  and  Allen  PB-X  (Baldwin-Llma-Hamilton)  cements. 
Although  both  cements  appeared  to  have  sufficient  bonding  strength  to 
pyrolytic  graphite  at  room  temperature,  they  failed  when  tested  at  high 
temperatures  up  to  1000*  F.  Since  this  bonding  problem  did  not  appear 
resolvable  within  a  reasonable  time  period,  efforts  in  this  direction 
were  discontinued.  Other  difficulties  were  also  anticipated,  such  as 
the  determination  of  the  corrections  necessary  for  the  thermal  expansion 
of  the  gage  itself  and  for  the  heating  of  the  lead  wires. 

Because  it  offered  certain  potential  advantages  over  the  electric- 
resistance  foil  gages,  the  Tuckerman  optical  strain  gage  system  was  finally 
selected  for  the  experimental  measurements. 

The  two  basic  units  of  the  system  are  the  extensometer  and  the 
autocollimator.  The  extensometer  consists  of  one  fixed  knife  edge  and 
one  rotating  knife  edge,  or  losenge,  on  which  there  is  a  mirrored  surface. 
The  two  knife  edges  are  placed  against  the  sample  under  test.  Deformation 
of  the  sample  causes  the  losenge  to  rock. 
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The  autocollimator  is  a  precision  telescope  consisting  of  a 
highly  corrected  objective  lens  system,  a  scale,  a  light  source  which 
projects  a  fiducial  spot  on  the  scale,  and  an  eyepiece. 

The  principle  by  which  the  two  units  work  is  illustrated  in 
the  schematic  diagram  of  Figure  5.  Light  emanating  from  the  autocollimator 
strikes  the  mirrored  surface  of  the  lozenge,  then  the  fixed  mirror,  and  is 
reflected  back  to  the  objective  lens  which  focuses  the  fiducial  spot  on 
the  scale.  The  position  of  this  spot  on  the  scale  is  a  function  of  the 
lozenge  rotation  caused  by  the  specimen  deformation.  Readings  obtained 
are  readily  converted  to  strain  values  by  a  simple  formula. 

The  extensometer  used  in  the  tests  was  a  standard  model  with 
a  1-inch  gage  length  and  a  0.4-inch  lozenge.  With  this  extensometer 
the  least  count  of  the  vernier  scale  readings  was  about  8  microinches 
per  inch. 

The  Tuckerman  optical  strain  gage  is  designed  primarily  for 

static  strain  measurements  at  steady  state  temperatures  and  is  recommended 
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for  use  up  to  a  temperature  of  500*  F.  The  limiting  condition  is 
excessive  heating  of  the  mirrored  surface  of  the  extensometer.  Because 
of  the  highly  transient  conditions  of  the  tests,  it  was  felt  that  the 
gage  might  be  usable  up  to  surface  temperatures  of  approximately  1000*  F 
without  damaging  the  mirrored  surface.  In  adapting  the  gage  for  measur¬ 
ing  the  transient  strains  in  a  pyrolytic  graphite  plate,  several  innova¬ 
tions  in  instrumentation  and  in  technique  were  necessary. 

A  specially  designed  stand  was  built  for  supporting  and  aligning 
the  autocollimator  above  the  test  section.  (See  Figure  2)  Alignment  of 
the  autocollimator  with  the  extensometer  presented  no  difficulties.  The 
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extensometer  was  placed  on  the  unheated  surface  of  the  plate  and  held  In 
place  by  a  spring-tensioned  wire  which  exerted  a  downward  force  of  about 
one  pound.  The  wire  was  easily  displaced  to  permit  immediate  removal  of 
the  extensometer.  The  unheated  surface  was  obviously  the  most  practical 
location  for  the  strain  measurements  since  the  lowest  temperatures  occurred 
there  and  it  was  readily  accessible  for  the  instrumentation. 

A  preliminaiy  test  on  a  l/8  inch  thick  plate  was  conducted  to  see 
if  the  response  of  the  Tuckerman  gage  could  be  followed  visually.  The 
results  were  favorable,  suggesting  that  a  sufficient  number  of  values  to 
define  a  curve  could  be  obtained  without  special  recording  equipment. 

Further  tests  were  conducted  to  measure  the  thermal  strains  in 
both  a  1/8  inch  and  a  1/4  inch  thick  plate.  A  push-button  switch  which 
transmitted  an  electrical  signal  to  a  recording  oscillograph  was  devised 
for  synchronizing  the  visual  strain  readings  with  a  time  scale.  Tempera¬ 
tures  at  the  unheated  surface  were  measured  near  the  leading  edge  by  a 
thermocouple  held  down  with  a  two  pound  spring  force.  The  length  of  time 
of  heating  was  gradually  increased  with  each  test  to  a  maximum  of  45 
seconds  for  the  l/8  inch  thick  plate  and  to  a  maximum  of  85  seconds  for 
the  1/4  inch  thick  plates.  Generally,  10  to  15  strain  readings  could  be 
taken  in  these  time  intervals.  A  maximum  temperature  of  1500*  F  was 
reached  on  the  unheated  surface  with  no  apparent  damage  to  the  mirrored 
surface  of  the  extensometer.  Excellent  reproduction  of  temperatures  from 
test  to  test  indicated  that  heating  rates  were  being  reproduced  quite  well. 

Considerable  scatter  in  the  strain  data  was  observed  in  the 
first  few  tests  of  each  plate.  Differences  were  also  noted  between  the 
strain  curves  of  the  "mandrel"  side  of  a  plate  and  the  "deposition"  side. 


However,  after  four  or  five  tests  on  each  side  of  the  plates,  consistent 
data  could  be  obtained.  This  behavior  was  attributed  to  the  presence  of 
residual  stresses  in  the  plates  which  apparently  were  eliminated  by  the 
annealing  effect  of  the  heating  and  cooling. 

During  the  course  of  these  preliminary  tests,  it  became  evident 
that  the  extensometer  was  being  subjected  to  considerable  heating.  It 
was  then  realized  that  the  resulting  thermal  expansion  of  the  extensometer 
was  registering  as  an  apparent  compressive  strain  and  seriously  Influencing 
the  strain  data. 

In  order  to  determine  the  correction  for  this  effect,  the 
thermal  expansion  of  the  extensometer  was  calibrated.  The  calibration 
was  performed  by  first  heating  the  extensometer  with  an  electric  blower- 
type  heater.  It  was  then  placed  on  an  unheated  pyrolytic  graphite  plate; 
and  as  it  cooled,  strain  readings  were  observed  through  the  autocollimator 
and  the  body  temperature  measured  by  a  thermocouple.  The  thermal  expansion 
of  the  pyrolytic  graphite  plate  was  believed  negligible  due  to  its  low 
thermal  expansion  coefficient  in  the  plane.  The  calibration  points  ob¬ 
tained  are  shown  in  Figure  6.  The  average  thermal  expansion  coefficient 
of  the  extensometer  was  found  to  be  7.30  microinches  per  inch  per  degree 
Fahrenheit. 

Of  course,  the  next  step  was  to  determine  the  amount  of  the 
temperature  rise  of  the  extensometer  during  a  test.  Measurements  showed 
that  this  temperature  rise  was  about  100°  F.  However,  this  temperature 
rise  was  not  uniform.  Temperature  differences  as  large  as  20*  F  were 
measured  between  locations  on  the  extensometer.  Since  there  was  no 
way  of  determining  the  "effective"  temperature  with  which  the  extensometer 
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expanded,  accurate  corrections  to  the  observed  strain  readings  were  not 
possible. 

A  solution  to  this  problem  was  seen  when  it  was  established  that 
about  807.  of  the  heat  transfer  to  the  extensometer  was  by  radiation.  The 
advantages  of  a  heat  shield  placed  between  the  extensometer  body  and  the 
plate  surface  were  then  immediately  obvious.  A  heat  shield  was  made  from 
thin  sheets  of  nichrome  separated  by  thin  insulating  sheets  of  mica. 

Holes  were  cut  in  the  shield  to  accommodate  the  knife  edge  and  lozenge. 
Wires  were  stretched  across  the  plate  to  hold  the  shield  in  position. 
Figure  7  shows  the  final  arrangement  of  the  instrumentation  in  the  test 
section. 

The  use  of  the  heat  shield  was  quite  effective  in  reducing  the 
radiation  heat  transfer.  The  temperature  rise  of  the  extensometer  was 
reduced  to  less  than  30*  F  and  the  maximum  temperature  difference  was 
reduced  to  about  3*  F,  corresponding  to  a  maximum  uncertainty  in  the 
correction  of  about  20  microinches  per  inch. 

3.5  Test  Procedure  for  Strain  Measurements 

The  uncertainties  in  the  observed  strain  readings  and  in  the 
corrections  to  these  readings  were  expected  to  cause  considerable  scatter 
in  the  data  points  for  a  particular  test.  With  this  in  mind,  a  number  of 
tests  were  conducted  on  both  the  l/8  inch  and  the  1/4  inch  thick  plates 
so  that  the  mean  trend  of  the  strains  for  each  plate  could  be  determined. 
The  strains  were  measured  along  the  center  line  of  each  plate  parallel 
to  the  flow  with  the  center  of  the  extensometer  located  approximately 
3*1/2  inches  from  the  trailing  edge.  A  few  additional  measurements  were 
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made  on  the  l/8  inch  thick  plate  with  the  center  of  the  extensometer 
located  approximately  1  inch  from  the  trailing  edge.  However,  theae 
latter  measurements  were  not  entirely  successful  and  will  be  discussed 
below.  Each  of  the  plates  tested  had  been  subjected  to  a  number  of 
previous  tests  so  that  any  residual  stresses  present  were  probably 
eliminated. 

The  procedure  followed  during  each  test  was  the  same.  The 
plate  was  placed  in  its  usual  position  (l/8M  above  and  1/2"  downstream 
of  the  nozzles)  in  the  test  section.  The  extensometer.  Instrumented 
with  two  thermocouples,  was  placed  on  the  plate  together  with  the  heat 
shield  and  held  in  place  by  the  spring-tensioned  wire.  It  was  necessary 
to  press  the  extensometer  against  the  plate  to  form  grooves  for  the 
knife  edges.  A  thermocouple  was  lowered  on  the  plate  to  measure  the 
temperature  near  the  leading  edge.  The  outputs  of  the  three  thermo¬ 
couples  were  recorded  on  the  Minneapolls-Honeywell  recording  oscillograph. 
Each  test  began  with  the  ignition  of  the  oxyacetylene  flame.  Care  had  to 
be  exercised  in  this  operation  because  an  explosive  ignition  would  usually 
knock  the  extensometer  off  scale  so  that  readings  could  not  be  taken.  The 
l/8  inch  thick  plates  were  heated  for  about  35  seconds  in  each  test  and 
the  1/4  inch  thick  plates  for  about  70  seconds.  The  visual  strain  readings 
observed  during  a  test  were  synchronized  with  a  time  scale  by  means  of 
the  push-button  switch.  At  the  termination  of  a  test,  the  extensometer 
was  imnedlately  removed  to  minimize  its  exposure  to  heating. 

During  the  early  part  of  a  test  the  strain  variation  was  so 
rapid  that  it  was  not  possible  to  read  the  vernier  scale.  However,  by 
interpolating  between  the  lines  of  the  primary  scale,  readings  could  be 
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taken  to  within  an  accuracy  of  about  ±  20  microinches  per  inch.  The 
strain  variation  decreased  considerably  during  the  latter  part  of  .the 
test  so  that  it  became  possible  to  read  the  vernier.  This  increase  in 
accuracy  was,  however,  offset  by  the  increase  in  the  uncertainty  of  the 
correction  for  the  extensometer  expansion  due  to  larger  temperature 
differences  in  the  extensometer. 

Considerable  difficulties  were  encountered  in  the  strain 
measurements  near  the  trailing  edge  of  the  l/8  inch  thick  plate.  The 
main  difficulty  was  due  to  flame  leakage  from  the  trailing  edge  of  the 
plate  which  considerably  increased  the  heat  transfer  to  the  extensometer. 
This  resulted  in  large  temperature  differences  in  the  extensometer  and 
invalidated  much  of  the  data  obtained.  Contact  between  the  heat  shield 
and  the  exhaust  duct  also  proved  to  be  troublesome.  The  vibrations  of 
the  duct  caused  by  the  flowing  gases  were  transmitted  through  the  heat 
shield  to  the  extensometer,  making  it  impossible  to  take  readings.  In 
spite  of  these  difficulties,  some  fairly  reasonable  data  were  obtained 


in  two  of  the  tests 
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4.0  EXPERIMENTAL  RESULTS 

A  large  amount  of  data  was  obtained  during  the  numerous  experimental 
runs.  All  of  the  data  were  analyzed  and  only  those  considered  to  be 
reasonably  valid  were  retained.  The  following  two  sections  present  the 
experimental  results  representing  the  best  data  obtained  from  the 
temperature  and  strain  measurements. 

4.1  Temperature  Distribution 

The  continuous  temperature  records  obtained  from  the  oscillograph 
for  both  the  l/8  inch  and  the  1/4  inch  thick  plates  are  shown  in  Figures 
8A  through  9E.  The  temperatures  are  plotted  as  functions  of  time  at 
various  depths  from  the  unheated  surface  and  each  figure  is  for  a  par* 
tlcular  distance  from  the  leading  edge.  These  temperature  records  were 
cross-plotted  in  two  ways.  First,  the  temperature  distributions  through 
the  thickness  were  plotted  at  various  times.  Figures  10  and  11  show 
these  distributions  for  both  plates  at  only  the  one  location  near  the 
trailing  edge  since  these  were  typical  and  were  the  only  ones  of  this 
type  required  for  the  calculation  of  the  thermal  stresses  and  strains. 

The  second  type  of  cross-plots  are  shown  in  Figures  12A  through  13G  from 
which  the  temperature  distributions  along  the  length  of  the  plates  can 
be  seen  at  various  depths.  Each  figure  is  for  a  particular  time  after 
heating.  Also  shown  in  these  figures  are  the  approximate  temperature 
curves  with  their  corresponding  equations  which  were  fitted  to  the 
experimental  data.  It  can  be  seen  that  these  fitted  curves  represent 
a  fairly  good  approximation  of  the  actual  temperature  distributions. 

Notice  that  the  temperatures  given  by  the  term  X„(X.)  are  just  the 
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temperatures  through  the  thickness  near  the  trailing  edge  which  were  shown 
in  Figures  10  and  11. 

Temperature  measurements  across  the  width  of  a  l/4  inch  thick 
plate  showed  that  the  maximum  variation  was  only  about  25*  F.^  There  was 
no  reason  for  believing  that  this  variation  was  very  much  larger  in  the 
l/8  inch  thick  plate.  Thus,  these  measurements  indicated  that  the  tempera¬ 
ture  distributions  in  the  plates  were  very  nearly  two-dimensional. 

4.2  Strain  Variation 

The  results  of  the  strain  measurements  taken  in  a  number  of  tests 
are  presented  in  Figures  14  and  15  where  the  strain  variation  with  time  is 
shown  for  both  the  l/8  inch  and  l/4  inch  thick  plates.  The  data  points 
which  are  plotted  have  been  corrected  for  the  thermal  expansion  of  the 
extensometer.  The  scatter  in  the  data  points  for  both  plates  is  approxi¬ 
mately  ±  20  microinches  per  inch,  except  for  the  data  points  of  tests  5 
and  6,  and  corresponds  to  the  experimental  error  (cf..  Section  3.5) 
expected  in  the  measurements.  Although  the  uncertainty  in  the  strain 
measurements  was  relatively  large,  it  is  seen  from  the  figures  that  the 
mean  trend  of  the  strain  variation  at  ■  3.5"  is  well  defined. 

The  Btralns  measured  near  the  trailing  edge  (x^  -  1.0")  of  the 
1/8  inch  thick  plate  are  represented  by  the  data  points  of  tests  5  and  6. 
Because  of  the  experimental  difficulties  previously  mentioned  (cf.,  Sec¬ 
tion  3.5),  these  strain  measurements  are  of  questionable  accuracy. 

However,  they  were  included  in  Figure  14  for  the  purpose  of  showing  the 
general  trend  of  strains  near  the  trailing  edge. 
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The  variation  of  the  strains  in  both  plates  is  seen  to  follow 
the  same  pattern.  In  each  case  the  initially  compressive  strains  reverse 
into  almost  linearly  increasing  tensile  strains.  This  behavior  is  com¬ 
pared  to  the  strain  variation  predicted  from  the  analytical  solution  in 
the  next  section. 
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5.0  DISCUSSION  OF  RESULTS 

It  is  axiomatic  in  experimental  stress  analysis  that  the  experimenter 
must  rely  directly  on  theoretical  analysis  for  his  knowledge  of  the  state 
of  stress  in  a  body.  Stresses  are  rarely  measured  directly;  whereas  strains 
are  almost  always  measured.  Therefore,  the  usual  test  of  the  validity  of 
a  theoretical  analysis  is  to  compare  the  experimentally  measured  strains 
with  the  strains  calculated  from  the  analytical  solution.  Then  if  agree¬ 
ment  is  found,  there  is  some  justification  for  believing  that  the  calcu¬ 
lated  stress  distribution  represents  the  actual  state  of  stress. 

In  Section  5.1  the  experimentally  measured  strains  are  compared  with 
the  strains  calculated  from  the  analytical  solution.  Then  in  Section  5.2 
the  calculated  stress  distributions  are  discussed.  The  effect  of  aniso¬ 
tropy  on  the  thermal  stresses  is  examined  in  Section  5.3.  Finally, 

Section  5.4  concludes  with  a  discussion  on  further  developments  suggested 
by  the  results  of  this  investigation. 

5.1  Comparison  of  Calculated  Results  with  Experimental  Data 

In  Figure  14  the  calculated  strain  variations  at  two  locations 
on  the  unheated  surface  of  the  1/8  inch  thick  pyrolytic  graphite  plate 
are  compared  with  the  experimental  data.  Reasonably  good  agreement 
(within  ±  20  micro  Inches  per  inch)  is  found  at  x^  -  3.5",  but  the 
comparison  is  inconclusive  at  x^  -  1"  due  to  the  limited  experimental 
data  available  at  this  location.  However,  the  predicted  strain  curve 
does  follow  the  general  trend  of  this  data.  Discrepancies  which  exist 
are  probably  due  to  the  uncertainty  in  the  strain  measurements. 

The  comparison  of  the  strain  variation  at  the  unheated  surface 
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of  the  l/4  inch  thick  pyrolytic  graphite  plate  with  the  experimental  data 
la  shown  in  Figure  15.  The  agreement  is  within  experimental  error  (1  20 
micro Inches  per  Inch)  up  to  about  50  seconds;  after  this  time  the  predicted 
strain  curve  gradually  diverges  from  the  data  to  a  maximum  difference  of 
80  microinches  per  inch  at  75  seconds.  No  definite  explanation  was  found 
for  this  deviation,  although  the  influence  of  the  constant  properties 
solution  (cf..  Section  2.2)  on  the  calculated  strains  is  suspected. 

In  regard  to  the  calculations,  it  should  be  mentioned  that  the 
calculated  strains  were  quite  sensitive  to  small  changes  in  the  assumed 
thermal  expansion  curve  in  the  ,,a"-direction  for  pyrolytic  graphite.  The 
thermal  expansion  curve  assumed  in  the  final  calculations  is  shown  in 
Figure  24.  The  use  of  this  curve  in  the  strain  calculations  was  found  to 
yield  the  best  fit  to  the  experimental  data  of  both  plates. 

5.2  Stress  Distributions 

The  behavior  of  the  strain  curves  presented  in  the  previous 
section  is  better  understood  by  first  examining  the  calculated  stress 
distributions  in  the  plates.  The  stress  distribution  of  will  not 
be  presented  here  since  it  is  not  very  different  from  the  distribution 
of  0^ 

Figures  16  and  17  show  typical  distributions  of  the  normal 
stress,  0^,  through  the  thickness  of  both  a  l/8  inch  and  1/4  inch 
thick  plate  at  various  locations  along  the  length.  The  stresses  at 
both  surfaces  of  a  plate  are  compressive,  the  larger  stress  occurring 
at  the  heated  surface.  The  compressive  stresses  near  the  surfaces  are 
exactly  balanced  by  the  tensile  stresses  deeper  in  the  plate.  The 
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stress  distribution  near  the  trailing  edge  (Xj/l  ■  l/6)  is  almost  solely 
due  to  the  temperature  distribution  T^Cx^.  Other  temperature  variations 
contribute  to  the  stress  distribution  at  distances  farther  along  the  length. 
These  are  seen  to  increase  the  stresses. 

Since  the  maximum  compressive  stresses  occur  at  the  heated  surface 
and  the  maximum  tensile  stresses  occur  near  the  midplane  (x^/h  -  0.5), 
the  stresses  at  both  of  these  locations  were  plotted  as  functions  of  time 
for  both  plates  in  Figures  18  and  19.  The  stresses  increase  during  the 
initial  heating  period,  but  gradually  approach  constant  values  as  the  tem¬ 
perature  gradients  in  the  plates  cease  changing.  This  behavior  is  in 
marked  contrast  to  the  strain  variations  in  the  plates  which  appear  to 
increase  without  limit. 

Figures  18  and  19  also  show  how  the  calculated  stresses  Increase 
with  increasing  distances  from  the  trailing  edge.  As  mentioned  above, 

x3)  are  responsible  for  this 
variation.  The  proportional  Increase  in  the  stresses  is  much  greater 
in  the  l/8  inch  thick  plate  than  in  the  l/4  inch  thick  plate.  This  is 
due  to  the  larger  lengthwise  temperature  gradients  which  occurred  in 
the  l/8  inch  thick  plates. 

Distributions  of  the  normal  stress,  Oy  perpendicular  to  the 
plane  of  the  plate  and  the  shearing  stress,  Oy  are  shown  in  Figures 
20  and  21.  These  are  typical  for  both  the  l/8  inch  and  l/4  inch  thick 
plates.  The  stresses  are  seen  to  be  of  negligible  magnitude.  Therefore, 
one  might  be  led  to  apply  the  usual  plane  stress  approximation  for  the 
two  cases  considered  here.  Upon  closer  examination,  however,  it  becomes 
clear  that  the  plane  stress  solution  would  lead  to  completely  erroneous 


temperature  variations  other  than  T^( 
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results  since  the  temperature  distribution  given  by  equation  (7)  does  not 
satisfy  simultaneously  both  the  assumptions  of  plane  stress  and  the  com¬ 
patibility  conditions.  This  is  more  easily  shown  by  inspection  of  either 
equations  (8)  or  equations  (31)  from  which  the  important  stresses  are 
seen  to  be  strong  functions  of  the  thickness  coordinate,  x^  (also  see 
Figures  16  and  17).  On  the  other  hand,  the  stresses  calculated  from 

ft 

plane  stress  are  assumed  independent  of  x^  ,  and  so  they  would  be 
totally  inconsistent  with  the  exact  results.  The  solution  may  be  slightly 
simplified  by  assuming  at  the  outset  that  o^,  c^,  and  a,.  are  zero, 

but  the  full  three-dimensional  equations  would  still  have  to  be  retained 
in  formulating  the  solution. 

The  strain  variation  in  the  plates  can  now  be  explained  on  the 
basis  of  the  information  presented  on  the  stress  distributions.  The 
strain,  e^,  at  the  unheated  surface  is  given  by  the  Hooke's  law 
relation  (4) 


« 


1 


13 


T 


where  the  stresses  and  the  thermal  expansion  are  evaluated  at  the  unheated 
surface.  Since  is  negligible  and  a 2  is  of  the  same  order  as  o^, 

the  strain  is  approximately 

(1  -  V12) 

*i  “  “*1^ —  a!+alT 


Note:  This  is  true,  at  least,  for  the  first-order  plane  stress  approxi¬ 
mation.  Higher-order  plane  stress  approximations  may  approach 

o 

the  exact  solution  for  certain  temperature  distributions  ,  but 
the  use  of  these  higher-order  approximations  probably  would  not 
simplify  the  solution. 
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Now  0^  at  the  unheated  surface  is  always  compressive  (See  Figures  16  and 

17),  and  at  low  temperatures  the  thermal  expansion,  Q^T,  is  negative. 

Thus,  it  is  seen  from  the  above  relation  that  the  strain  will  be  compressive 

during  the  early  heating  period  (See  Figures  14  and  15).  At  later  times 

0^  becomes  nearly  constant  (Figures  18  and  19),  while  O^T  becomes 

positive  and  increases  almost  linearly  with  the  temperature  rise.  The 

thermal  expansion  eventually  overtakes  the  compressive  strains  produced 

by  0^  and  then  dominates  the  strain  variation.  The  tensile  strain 

variation  is  due  to  this  effect. 

The  stress  levels  produced  in  the  pyrolytic  graphite  plates 

are  relatively  low  compared  to  the  ultimate  tensile  strength  of  about 
20 

20,000  psl  measured  for  pyrolytic  graphite.  Figures  18  and  19  show 

that  maximum  stresses  of  -1400  psl  and  -3000  psl  were  produced  in  the 

1/8  inch  and  l/4  inch  thick  plates,  respectively.  However,  one  should 

bear  in  mind  that  the  plates  were  subjected  to  a  relatively  low  heating 

rate  of  about  40  Btu/sec-ft  ,  and  that  higher  heating  rates  would  lead 

to  more  severe  temperature  gradients,  which  in  turn  would  increase  the 

thermal  stresses.  When  one  considers  the  fact  that  a  heating  rate  of 
2 

1000  Btu/sec-ft  is  not  uncommon  during  the  re-entry  of  space  vehicles, 
the  increased  significance  of  thermal  stresses  is  appreciated. 

5.3  Effect  of  Anisotropy 

The  effect  of  the  degree  of  anisotropy  on  the  thermal  stresses 
in  a  transversely  anisotropic  plate  can  be  conveniently  examined  by  con¬ 
sidering  the  particular  temperature  distribution  given  by 


T  -  bx* 
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From  equation  (31)  the  normal  stress  arising  from  this  temperature 

distribution  la 


°1  ’ 


O^Ej^bh 


1  -  V 


12 


En  a 

v13  E  +  a 
13  e33  a1 


The  maximum  value  of  the  terms  in  the  brackets  is  -  0.1667.  The  maximum 
stress  can  therefore  be  written  in  the  following  form: 


<1 


v12) 


a1E11bh 


-  0.1667 


13  E33  +  / 


This  equation  shows  that  the  effect  of  elastic  anisotropy  is  contained  in 

the  ratios  ^n/^33  and  a^/cXy  In  order  to  examine  this  effect  by  itself, 

2 

it  is  convenient  to  specify  that  the  ratio  (1  -  Vj^A-jE*  p  the  tempera- 

2 

ture  distribution  T  «  bx^  ,  and  the  plate  thickness  h  all  remain 
constant.  Then  the  maximum  stress  will  simply  increase  linearly  with 

r  eh  a3 1 

the  anisotropic  parameter  |^V13  g — +  —  J  as  shown  in  Figure  22.  Values 
of  this  parameter  corresponding  to  a  pyrolytic  graphite  plate  and  to  an 
isotropic  plate  are  indicated  on  the  figure.  The  larger  value  for  the 
pyrolytic  graphite  plate  is  due  to  the  large  value  of  a3  compared  to 
0^.  It  is  interesting  to  note  that  the  maximum  stress  in  the  pyrolytic 
graphite  plate  is  more  than  10  times  the  maximum  stress  in  an  isotropic 
plate. 

The  foregoing  analysis  may  seem  somewhat  superficial  since,  in 
general,  the  thermal  properties  would  vary  with  changes  in  the  elastic 
properties,  resulting  in  different  temperature  distributions  for  a  given 
heating  rate.  For  example,  instead  of  comparing  the  pyrolytic  graphite 
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plate  to  the  Isotropic  plate  with  identical  thermal  properties,  it  might 

have  been  compared  to  an  isotropic  plate  of  ordinary  graphite.  In  this 

2 

case,  the  temperature  gradients  of  the  type  T  ■  bx^  would  be  much  smaller 
in  the  pyrolytic  graphite  plate  than  in  the  ordinary  graphite  plate  due  to 
its  much  higher  thermal  conductivity  in  the  plane.  The  higher  stresses  in 
the  pyrolytic  graphite  plate  caused  by  elastic  anisotropy  might  then  be  just 
compensated  by  the  smaller  temperature  gradients,  so  that  the  maximum 
stress  in  the  pyrolytic  graphite  plate  may  be  about  the  same  as  that  in  the 
ordinary  graphite  plate.  Thus,  if  one  were  using  the  comparison  of  thermal 
stresses  as  the  criterion  for  selecting  a  material  for  application,  it 
would  be  important  to  account  for  the  difference  in  temperature  distribu¬ 
tions.  However,  in  practice,  the  heat  shielding  capability  of  a  material 
is  usually  the  governing  criterion,  in  which  case  there  would  be  no  ques¬ 
tion  about  the  superiority  of  pyrolytic  graphite  over  ordinary  graphite. 
Moreover,  the  purpose  of  the  above  analysis  was  simply  to  provide  basic 
insight  into  the  mechanism  of  elastic  anisotropy.  Although  it  may  be 
more  realistic  physically  to  do  so,  the  inclusion  of  temperature  effects 
would  merely  tend  to  obscure  the  analysis  because  a  general  relationship 
does  not  appear  likely  to  exist  between  the  elastic  properties  of  materials 
and  their  thermal  properties. 

5.4  Recommenda  t ions 

In  this  section  a  few  suggestions  are  made  for  extending  the 
theoretical  analysis,  and  some  recommendations  are  forwarded  for  improv¬ 
ing  the  experimental  technique. 

The  generalization  of  the  solution  presented  in  Section  2.3  to 
any  two-dimensional  temperature  distribution,  T  ■  TCx^x-j),  does  not 
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appear  to  present  any  basic  difficulties.  It  should  be  noted  that  the 
derivation  of  equation  (16),  the  basic  governing  differential  equation, 
was  not  based  on  the  assumption  of  a  particular  two-dimensional  tempera¬ 
ture  distribution.  In  principle,  a  solution  of  (16)  could  be  derived 
for  any  two-dimensional  temperature  distribution  by  repeated  integration 
of  equations  (19).  However,  solutions  of  the  compatibility  equations, 
(6d)  and  (6f),  for  the  stresses,  o^  and  0g,  may  present  additional 
difficulties. 

There  also  appears  to  be  a  good  possibility  of  extending  the 
solution  to  an  orthotropic  plate.  Basically,  this  would  simply  Involve 
some  changes  in  the  property  coefficients  of  the  transformed  compati¬ 
bility  conditions  (6).  Such  changes  probably  would  not  introduce  any 
serious  difficulties,  but  the  problem  has  not  been  investigated. 

Although  the  500*  F  extensometer  (cf..  Section  3.4)  proved 
to  be  adequate  for  the  purposes  of  this  investigation,  it  may  leave  much 
to  be  desired  for  strain  measurements  at  higher  temperatures  and  under 
highly  transient  conditions.  Considerable  improvement  in  this  technique 
could  probably  be  achieved  by  the  development  of  an  extensometer  designed 
specifically  for  such  conditions.  This  extensometer  should  incorporate 
the  following  features: 

1.  Knife  edges  constructed  from  a  material  of  low  thermal 
conductivity  which  can  withstand  very  high  temperatures. 

2.  Body  constructed  from  a  material  of  high  thermal  con¬ 
ductivity  and  low  thermal  expansion. 

3.  A  radiation  heat  shield  integral  with  the  body  of  the 


extensometer 
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These  features  would  reduce  the  heat  transfer  to  the  body  of  the  extenso- 
meter  and  also  reduce  the  uncertainties  Involved  in  the  correction  for  the 
thermal  expansion  of  the  gage  Itself.  Although  an  extensometer  designed 
for  use  at  temperatures  up  to  1000*  F  has  beer  described  in  Reference  23, 
it  is  evidently  not  yet  available  commercially. 

In  addition  to  these  modifications  of  the  extensometer,  a  means 
of  recording  the  strain  readings  by  other  than  visual  observations  would 


be  desirable 
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6.0  CONCLUSIONS 

The  results  of  this  investigation  lead  to  the  following  conclusions: 

1.  The  general  equations  of  thermoelasticity  were  developed  for  a 
transversely  anisotropic  material  such  as  pyrolytic  graphite. 

The  equations  were  solved  exactly  for  a  particular  two-dimensional 
temperature  distribution  in  a  flat  plate  to  yield  a  solution  for 
the  stress-strain  distribution. 

2.  A  technique  was  described  for  measuring  transient  strains  up  to 
a  temperature  of  1500*  F  using  a  Tuckerman  optical  strain  gage. 

3.  The  strains  measured  in  heated  pyrolytic  graphite  plates  by  this 
technique  were  generally  in  agreement  within  experimental  error 
(1  20  microinches  per  inch)  with  the  strains  calculated  from  the 
analytical  solution. 

A.  The  calculated  stress  distributions  showed  that  only  the  stresses 
(0^  and  02)  in  a  direction  parallel  to  the  plane  of  a  plate  were 
of  significant  magnitude.  Although  the  shearing  stresses  (0^, 

0y  and  0g)  and  the  stress  normal  to  the  plane  (o^)  were  negli¬ 
gible,  the  plane  stress  approximation  was  found  to  be  inapplicable. 
5.  The  thermal  stresses  in  a  transversely  anisotropic  plate  were 
shown  to  increase  linearly  with  the  anisotropic  parameter, 

r  En  a3i  2 

LV13  E —  +  cT”J  *  *or  temPerature  distribution  T  •  bx^  . 

For  this  case,  the  maximum  stress  in  a  pyrolytic  graphite  plate 
was  found  to  be  more  than  10  times  that  in  an  isotropic  plate. 
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SECTION  A-A 


FIG  3  DETAIL  OF  TEST  SECTION  ASSEMBLY 


HYD  7669 


57 


FIGURE  4  A  -  SCHEMATIC  OF  THERMOCOUPLE  HOLDING  DEVICE  AND  CIRCUIT 


FIGURE  4  b  -  TEMPERATURE  MEASUREMENT  SYSTEM 
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FIG.  6  CALIBRATION  OF  TUCKERMAN  EXTENSOMETER  (0.4  IN  LOZENGE  ) 
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FIG.  7  STRAIN  MEASUREMENT  SYSTEM 
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FIG. 88  TEMPERATURE  RISE  AT  VARIOUS  DEPTHS 

IN  A  '/*  INCH  THICK  PYROLYTIC  GRAPHITE 
PLATE  I  ^  INCHES  FROM  LEADING  EDGE 
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FIG.  9C  TEMPERATURE  RISE  AT  VARIOUS  DEPTHS 

IN  A  '/♦  INCH  THICK  PYROLYTIC  GRAPHITE 
PLATE  \%  INCHES  FROM  LEADING  EDGE 
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FIG.  10  TEMPERATURE  DISTRIBUTION  AT  VARIOUS  TIMES 
THROUGH  I/*  INCH  THICK  PYROLYTIC  GRAPHITE 
PLATE.  5  INCHES  FROM  LEADING  EDGE 
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DISTANCE  FROM  TRAILING  EDGE  OF  PLATE- INCHES 


FIG  .  12  A  TEMPERATURE  DISTRIBUTION  ALONG  LENGTH  OF 

'/e  INCH  thick  pyrolytic  graphite  plate  at 
VARIOUS  DEPTHS  AFTER  5  SECONDS  OF  HEATING 
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FIG.I2B  TEMPERATURE  DISTRIBUTION  ALONG  LENGTH  OF 
INCH  THICK  PYROLYTIC  GRAPHITE  PLATE  AT 
VARIOUS  DEPTHS  AFTER  10  SECONDS  OF  HEATING 
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FIG.I2C  TEMPERATURE  DISTRIBUTION  ALONG  LENGTH  OF 
INCH  THICK  PYROLYTIC  GRAPHITE  PLATE  AT 
VARIOUS  DEPTHS  AFTER  15  SECONDS  OF  HEATING 
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FIG.  12  D  TEMPERATURE  DISTRIBUTION  ALONG  LENGTH  OF 
'/a  INCH  THICK  PYROLYTIC  GRAPHITE  PLATE  AT 
VARIOUS  DEPTHS  AFTER  20  SECONDS  OF  HEATING 
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FIG  .I2E  TEMPERATURE  DISTRIBUTION  ALONG  LENGTH  OF 
V%  INCH  THICK  PYROLYTIC  GRAPHITE  PLATE  AT 
VARIOUS  DEPTHS  AFTER  30  SECONDS  OF  HEATING 
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FIG.  I3A  TEMPERATURE  DISTRIBUTION  ALONG  LENGTH  OF 
‘/♦INCH  THICK  PYROLYTIC  GRAPHITE  PLATE  AT 
VARIOUS  DEPTHS  AFTER  10  SECONDS  OF  HEATING 
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FIG.  13  B  TEMPERATURE  DISTRI BUTION  ALONG  LENGTH  OF 
'/«•  INCH  THICK  PYROLYTIC  GRAPHITE  PLATE  AT 
VARIOUS  DEPTHS  AFTER  20  SECONDS  OF  HEATING 
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FIG.  13 C  TEMPERATURE  DISTRIBUTION  ALONG  LENGTH  OF 
'/4  INCH  THICK  PYROLYTIC  GRAPHITE  PLATE  AT 
VARIOUS  DEPTHS  AFTER  30  SECONDS  OF  HEATING 
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FIG  1 3D  TEMPERATURE  DISTRIBUTION  ALONG  LENGTH  OF 
'A  INCH  THICK  PYROLYTIC  GRAPHITE  PLATE  AT 
VARIOUS  DEPTHS  AFTER  40  SECONDS  OF  HEATING 
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FIG.  13  E  TEMPERATURE  DISTRIBUTION  ALONG  LENGTH  OF 
%  INCH  THICK  PYROLYTIC  GRAPHITE  PLATE  AT 
VARIOUS  DEPTHS  AFTER  50  SECONDS  OF  HEATING 
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FIG  13 F  TEMPERATURE  DISTRIBUTION  ALONG  LENGTH  OF 
■/♦  INCH  THICK  PYROLYTIC  GRAPHITE  PLATE  AT 
VARIOUS  DEPTHS  AFTER  60  SECONDS  OF  HEATING 
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FIG.  13  G  TEMPERATURE  DISTRIBUTION  ALONG  LENGTH  OF 
'/♦  INCH  THICK  PYROLYTIC  GRAPHITE  PLATE  AT 
VARIOUS  DEPTHS  AFTER  70  SECONDS  OF  HEATING 
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FIG.  14  COMPARISON  OF  EXPERIMENTAL  DATA  WITH  CALCULATED 
STRAINS  AT  UNHEATED  SURFACE  OF  1/8  INCH 
THICK  PYROLYTIC  GRAPHITE  PLATE 
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FIG.  15  COMPARISON  OF  EXPERIMENTAL  DATA  WITH  CALCULATED 
STRAINS  AT  UNHEATED  SURFACE  OF  1/4  INCH 
THICK  PYROLYTIC  GRAPHITE  PLATE 
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FIG.  16  DISTRIBUTION  OF  <7i  THROUGH  '/SINCH  THICK  PLATE 
AFTER  20  SECONDS  OF  HEATING 
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FIG.  17  DISTRIBUTION  OF  C7]  THROUGH  */4  INCH  THICK  PLATE 
AFTER  40  SECONDS  OF  HEATING 


FIG.  18  VARIATION  OF  (77  AT  HEATED  SURFACE  AND 
MIDPLANE  OF  1/8  INCH  THICK  PLATE 
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FIG.  20  DISTRIBUTION  OF  0~^  THROUGH  '/4  INCH  THICK 
PLATE  AFTER  40  SECONDS  OF  HEATING 
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FIG. 21  DISTRIBUTION  OF  THROUGH  1/4  INCH  THICK 
PLATE  AFTER  40  SECONDS  OF  HEATING 


FIG.  22  VARIATION  OF  THERMAL  STRESS  WITH  DEGREE  OF  ANISOTROPY 
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FIG.  23  MODULUS  OF  ELASTICITY  OF  PYROLYTIC  GRAPHITE 
IN  THE  "  a"  DIRECTION  v*  TEMPERATURE 


THERMAL  EXPANSION,  <*,T  =  *L/L-  xIO4  IN/IN 
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FIG.  24  THERMAL  EXPANSION  OF  PYROLITIC 
GRAPHITE  IN  THE  "  &"  DIRECTION 
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THERMAL  EXPANSION, o<  ,T  =  Al/L  -  x  I03  IN /IN 
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FIG.  25  THERMAL  EXPANSION  OF  PYROLYTIC 
GRAPHITE  IN  THE  "C  "  DIRECTION 
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REDUCTION  OF  INDEPENDENT  ELASTIC  CONSTANTS  FOR  VARIOUS  TYPES  OP  SYMMETRY 

For  the  general  anisotropic  solid,  the  generalised  Hooke's  law  Includ¬ 
ing  thermal  strain  Is  expressed  by** 


*11  "  3*^  "  allffll  +  a12*22  +  *13*33  +  a14°23  +  a15*13  +  *16*12  +  °1T 

*22  "  lx^  "  a12ffll  +  a22*22  +  *23*33  +  *24*23  +  a25°13  +  *26*12  +  °2T 

*33  "  3x^  “  *13*11  +  a23°22  +  *33*33  +  *34*23  +  a35°13  +  *36*12  +  a3T 

*23  “  +  “  a14fl,ll+a24°22+a34033+a44023+a45013+a46012  +  a4T 

*13  "  +  "  *15*ll+*25*22+*35*33+*45*23+a55*13+*56*12  +  a5T 

*12  “  ^  +  ISj  "  a16°ll+a26022+a36033+a46023+a56013+a66°12  +  V 


>(A.D 


It  Is  convenient  to  Introduce  the  following  notation: 


*11  “  *1 


*11  "  °1 


*22  “  *2 


°22  *  °2 


*33  “  *3 


*33  "  *3 


*23  “  *4 


*23  "  *4 


*13  "  *5 


*13  “  *5 


*12  "  *6 


*12  "  *6 


Then  equations  (A.l)  may  be  abridged  as: 

*1  "  alk*k  +  aiT  (1  -  1,  2 . 6) 
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The  Isothermal  elastic  coefficients,  aik,  are  defined  as 


where  the  subscript  T  denotes  constancy  of  temperature  and  subscript 
{ 0}  denotes  constancy  of  all  0^  other  than  0^.  Since  the  indices 
of  a^  can  take  six  values,  it  comprises  a  6  x  6  matrix.  However, 
the  symmetry  condition,  a^  -  a^,  reduces  the  number  of  independent 
coefficients  to  twenty-one.  Then  the  general  matrix  can  be  written  as 


all 

*12 

a13 

*14 

*15 

*16 

•  •  e 

a22 

a23 

*24 

a25 

*26 

•  •  e 

see 

a33 

*34 

*35 

*36 

•  #  • 

see 

see 

*44 

*45 

*46 

•  •  a 

see 

see 

see 

*55 

*56 

#  •  e 

see 

see 

see 

see 

*66 

The  are  symmetric  about  the  diagonal. 

12 

The  coefficients  of  thermal  strain,  are  defined  by 

ai  "  [^){o}  (i"  1#  6> 

where  the  subscript  (0}  denotes  constancy  of  stresses. 

The  nusiber  of  Independent  elastic  constants,  a^»  of  a  solid  reduce 
considerably  when  the  solid  exhibits  elastic  symmetry  with  respect  to  a 
plane  or  about  an  axis.  A  solid  is  said  to  exhibit  elastic  symmetry 
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I 


I 


with  respect  to  a  transformation  if  the  elastic  constants  are  Invariant 
in  this  transformation.  Consider  the  transformation  of  one  system  of 
axes  Xj,  x2,  x^  to  another  system  of  axes  xj,  x£,  x^  .  Suppose  this 
transformation  is  connected  by  the  orthogonal  scheme 


X1 

X2 

x3 

X1 

X1 

J1 

“l 

nl 

X* 

2 

J2 

«2 

n2 

X  * 

3 

J3 

m3 

n3 

where  J,  m,  n  denote  direction  cosines.  It  can  be  shown  that  the  corres¬ 
ponding  transformation  of  stresses  and  strains  are  given  by** 

2  2  2 

o'  -  i‘o1  +  m^02  +  n‘cr3  +  ZiDjUjO^  +  20^ 


°6  “  ^1^2ffl  +  “i0^^  *  nln2°3  +  ^mlnl  +  m2ni^  a4 

+  (nji2  +  ngjj)  o5  +  (JjBVj  +  J^) 


*1  "  il*l  +  “l*2  +  nl*3  +  mlnlf4  +  nlil*5  +  ^lml*6 


*4  "  2J2J3*1  +  2o2B3€2  +  2n2n3*3  +  (n2m3  +  m3n2)  *4 

^  (njij  ^  n3i2)  Cj  +  (^2m3  ^3^*2^ 

eeeee  e  e  •  e  • 

As  a  specific  example,  consider  the  case  of  one-plane  symetry  with 
respect  to  the  x^-plane  such  that 


101 


*1  -  x’  ,  x2  -  x'  ,  x3  -  -  *3 


The  direction  cosines  of  this  transformation  are  as  follows 


-1 


Transformations  of  stresses  and  strains  give 

°1  “  ®i  *  *[  "  *i  <*  “  2»  3»  6> 

°4  ~«4  •  «4  “  ”*4  •  °5  "  -°5  *  *5  "  ’*5 

For  the  x|,x£,x^  system,  the  Hooke's  law  relation  for  the  strain  Is 

‘i  -  ‘a0!  +  ‘a'i  +  ‘a0}  +  *i*°«  +  ‘is's  +  ‘le's  +  °iT 


Transforming  to  the  x^x^x^  system,  the  strain  becosMs 


*1  “  all°l  +  a12°2  +  *13°3  '  a14°4  "  *15a5  +  a16°6  +  °1T 

However,  since  the  a^  are  invariant,  the  strain  can  be  written  without 
reference  to  the  xj^x^x^  system  as 


*1  "  allal  +  a12a2  +  a13®3  +  a14C4  +  a15°5  +  a16°6  +  “lT 


Comparison  of  the  two  equations  yields 


a14  "  *15  "  0 
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Examining  <2*  t^,—  in  turn,  the  following  oust  also  hold: 


l24  “  *25  "  *34  "  *35  “  *46  “  *56  "  0 


Thus  the  elastic  matrix  for  one -plane  symmetry  with  respect  to  the 
plane  can  be  written  as 


*n 

*12 

*13  0 

0 

*16 

*22 

*23  0 

0 

*26 

*33  0 

0 

*36 

*ik  " 

*45 

*46 

0 

*55 

0 

*66 

By  following  the  same 

procedures,  it  can 

be  shown 

that 

(orthotropic)  symmetry 

the  matrix  reduces 

to 

*n 

*12 

*13  0 

0 

0 

*22 

*23  0 

0 

0 

*33  0 

0 

0 

■ 

Ik 

*44 

0 

0 

*55 

0 

13  Independent 
constanta 
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9  independent 
constants 


And  for  the  solid  which  is  transversely  isotropic  in  the  XjXj^plane 
such  that  the  x^-axis  is  an  axis  of  syonetry,  the  matrix  reduces  to 


t 
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*11 

*12 

*13 

0 

0 

0 

*11 

*13 

0 

0 

0 

*13 

0 

0 

0 

*44 

0 

0 

*66 


5  Independent 
constants 
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CALCULATIONS 

i 

II.A  Properties  of  Pyrolytic  Graphite 

Although  extensive  work  has  been  done  on  evaluating  the  properties 

1 

of  pyrolytic  graphite,  much  of  the  available  data  is  only  preliminary  and 

some  of  the  properties,  such  as  the  shearing  modulus  g44.  have  not  been 

measured.  In  some  cases  the  properties  can  only  be  estimated;  while  in 

other  cases  a  reasonable  choice  must  be  made  between  conflicting  data. 

The  modulus  of  elasticity  in  the  "a"  direction  is  plotted  as  a 

function  of  temperature  in  Figure  23.  These  data  were  taken  from  Reference 

22 

,  21.  More  recent  data  indicate  that  the  "a"  direction  modulus  may  be 

as  large  as  10  x  10**  psl  in  compression.  It  is  also  indicated  that  the 
"c"  direction  modulus  may  exhibit  the  same  type  of  behavior.  The  "c" 
direction  modulus  in  tension  has  been  measured  to  be  1.5  x  106  psl,  and 
in  compression  it  has  been  measured  to  be  16  x  10**  psi. 

The  thermal  expansion  in  the  "a"  direction  and  in  the  "c" 
direction  is  plotted  against  temperature  in  Figures  24  and  25,  respec¬ 
tively.  These  data  were  also  taken  from  Reference  21.  The  difference 
in  the  thermal  expansions  of  the  top  and  bottom  layers  is  attributed 
to  the  presence  of  residual  stress.  The  solid  curve  in  Figure  24  is 

the  average  curve  used  in  the  numerical  calculations. 

21 

A  value  of  0.24  was  reported  to  have  been  measured  for  the 

Poisson's  ratio  in  the  basal  plane  (v^)*  This  ratio  in  the  lateral 

direction  (V..)  has  not  been  measured  yet. 

>  J 

The  value  of  the  shear  modulus  G44  is  unknown,  but  judging  from 
the  laminate-type  structure  of  pyrolytic  graphite,  it  seems  reasonable 


to  assume  that  greater  rigidity  will  be  exhibited  to  shearing  forces 
normal  to  the  basal  plane  than  to  shearing  forces  parallel  to  this 
plane.  Then  the  following  inequality  is  noted: 
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^  11 
G44  >  G11  "  2(1  +  V12) 

For  calculation  purposes,  the  following  approximation  was  made: 

G44  “  2G11  "  (1  +^2) 

II. B  Calculation  of  Stresses  and  Strains 

The  stresses  and  strains  are  given  by  equations  (8)  and  (9)  which 
account  for  property  variations,  and  by  equations  (31)  and  (32)  which 
are  based  on  constant  properties.  Numerical  calculations  were  carried 
out  for  the  strain,  C^,  at  the  unheated  surface  (Z»  0)  and  for 
the  stresses,  9^,  o^,  and  Og. 

In  evaluating  equations  (8)  and  (9)  for  and  9^  (e..  and 

9^  are  zero)  the  elastic  modulus  was  taken  as  a  constant,  along  with 
the  Poisson's  ratio,  since  it  does  not  vary  much  from  room  temperature 
to  3000*  F  (cf.,  Figure  23).  The  following  values  were  assumed: 

Eh  »  5.0  x  106  psi 

Vl2  -  0.24 

The  Integrals  and  were  evaluated  numerically  by  dividing  the 
thickness  of  each  plate  into  20  equal  increments,  The  nusierlcal 
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t 


i t 


i 


express Ions  for  end  o^  eppeer  as  follows: 


1  20  -  3  20  -  1 

«l<x3-°)  -  20  L  0^  T  -  10h  L  T  (*3  -  2  h> 

1-1  i-1  i 


(in/ In) 


°1<X3>  "  ^i°-*0^24  (  ’  alT(x3>  +  20  ai  * 


3(x,  -  r  h)  20  _  .  I 

+  — *~2  —  L  aiT  <*3  - 1  h>  r  <P»*> 

5Yi  1-1  1  z  J 


where: 

x.  •  average  distance  of  increment  from  x1 -axis 
J1  1 

Ot.T  -  thermal  expansion  at  average  temperature  T  of 

A  th 

1  inclement 

T(x3)  ■  temperature  at  x3> 


The  temperature  curves  used  in  these  calculations  are  shown  in  Figures 
10  and  11.  The  thermal  expansion,  a^T  ,  was  obtained  from  the  solid 
curve  of  Figure  24. 

In  evaluating  equations  (31)  and  (32)  for  o^,  Oy  and  Oy 
the  following  average  values  of  the  properties  were  assumed: 


:11 

- 

5.0  x  10^  psl 

'33 

■1 

15.0  x  106  psl 

12 

■ 

0.24 

13 

- 

0.30 

'44 

- 

4.0  x  10^  psl 

» 
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1.00  X  10"6  *F_1  for  the  l/8"  thick  plete 

a,  -  ,  , 

0.75  x  10  #F  for  the  l/4"  thick  plate 

a3  -  13.0  x  10"6  #f"1 


The  larger  compressive  value  of  E^3  was  chosen  because  o^  is  always 
in  compression  (See  Figure  20).  The  larger  value  of  oc^  was  assumed 
for  the  l/8  inch  thick  plate  because  of  the  slightly  higher  temperatures. 

The  parameters  defined  by  equation  (17)  were  then  calculated  to 
have  the  following  values: 


P  -  1.063 

f  6.58  pai/*F 

H  -  < 

L  4.93  psi/*F 

69.5  psi/*F 

C  -  « 

69.4  psi/*F 


1/8"  plate 
1/4"  plate 
1/8"  plate 
1/4"  plate 


7  was  not  calculated  since  it  does  not  appear  in  the  final  equations. 

Substitution  of  these  parameter  values  in  equations  (31)  and  (32) 
yielded  the  following  numerical  expressions  for  o^,  o^,  and  o^: 

For  the  1/8"  thick  plate: 


61| 


x  3-0 
xj-l" 


-  [2.245  c  x  10‘2  +  1240  b  +  a]  x  10’6  in/in 


t.i  -  «.i  +  [-  5.50  c  x  10'2  +  10.35  b  +  2.5  a]  x  10-6  (in/in) 

A|x,-0  i'x-“0 

xJ-3.5"  xj-l" 
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ax  -  3.70  CX2  (  -  £+  Z-Z2)  +  1.086  b  (  -  i  +  z  -  Z2  ) 

+  1.696  x  10"2  c  (  -  ^i  +  yfz-lz4)*  negligible  (psi) 

tern 

For  the  l/4"  thick  plate: 

e.i  -  [-  10.91  c  x  10'2  +  2.63  a]  x  lO-6  (ln/ln) 

1 ' x_»0 
Xj-3.5" 

Ox  m  11.09  c  X2  (  -  £  +  Z  -  Z2  )+  4.34  b  (  -  |  +  Z  -  Z2  ) 

+  0.271  c  (  “35  +  xf  ^  )+  negligible  (pal) 

term 

”3Uo..c'4-34(-£z2  +  5z3-K)  <-l> 

a. i  -  -  17.3  x(  -  4  Z+  Z2  -  4  Z3  I  (pal) 

3 « T«40»ec  »  J  J  ' 


The  constanta  a,  b,  and  c  are  the  coefficients  of  the  temperature 
curves  in  Figures  12. A  through  13.6.  The  total  strain  and  the  total 
stresses,  o^,  o^,  and  0^  wars  found  by  adding  the  results  of  the 

separate  calculations. 
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